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Total and partial cross sections for breakup of ground rovibronic state of by photon impact 
are calculated using the exact nonadiabatic nonrelativistic Hamiltonian without approximation. 
The converged results span six orders of magnitude. The breakup cross section is divided into 
dissociative excitation and dissociative ionization. The dissociative excitation channels are divided 
into contributions from principal quantum numbers 1 through 4. For dissociative ionization the 
kinetic energy sharing is calculated using a formally exact expression. These results are compared 
with approximate expressions, and it is shown that the Born-Oppenheimer result is surprisingly 
accurate, whereas using Born-Oppenheimer final states to extract the cross sections from the full 
nonadiabatic wave function produces pathologies near threshold. 

PACS numbers: 31.15.-p, 33.80.Eh, 31.15.xv 



The H^ cation is the smallest molecule, and one that 
is relevant in contexts ranging from interstellar chem- 
istry [J to fusion reactors, and as such is well studied 
in the literature. It provides the one-electron archetype 
for fundamental processes such as dissociative recombi- 
nation Due to its size, it is tractable to include 
nonadiabatic effects in calculating its dynamics [H, H|. 
Such studies provide insight into how coupled electronic 
and nuclear dynamics may be manipulated by laser light, 
and in recent years the interest in strong field and ultra- 
fast physics has led to many experimental and numeri- 
cal studies on this topic including accurate descriptions 
of charge localization [|| El , various processes in strong 
fields [7 I , and others [12l . [l3| . Benchmark calcula- 
tions of bound state rovibronic energies have been given 
in Refs. pJ4ll. 

Given the degree to which it is studied, it is surprising 
that no ab initio calculations of its one-photon, Fermi's 
golden rule breakup cross section that treat the nuclear 
and electronic degrees of freedom on the same footing 
have been published. The process represents one of the 
three fundamental coulomb breakup problems, the others 
including double ionization of helium, a complete calcula- 
tion of which was reported in 2005 [2(| ■ The present cal- 
culations are less computationally demanding than those 
on double ionization, but are more difficult to analyse. 



HAMILTONIAN AND BASIS 



conditions. 

The exact nonrelativistic Hamiltonian may be writ- 
ten m 



H = — 



1 



1 



1 



1 



1 



V + „ + — x 



T, 



2p e R 2 R ta tb Hr 

J (J + 1) - 2J 2 + J+f- + J-1+ + P 



2R 2 



2 x 1836.152701 
2 x 1836.152701+1 



1 



fJ-R 



-1836.152701 



(1) 



with the interparticle distances R, ta-, and rs, with \/ 2 
the laplacian in the electronic coordinates, and J z the 
projection of angular momentum (total, J, and elec- 
tronic, I; the projection of nuclear angular momentum 
is zero) upon the bond axis, conjugate to the third Eu- 
ler angle 7. Except for \/ 2 operators are denoted with 
hats and scalars have no hats in this equation. I calcu- 
late i? 5 / 2 times the wave function, such that the nuclear 
kinetic energy may be written 
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in which expression 
It may be shown that 
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The calculations employ an implementation [2lM24l | of 
the full nonadiabatic Hamiltonian in prolate spheroidal 
coordinates [J, l23j ■ The nuclear basis set is identical to 
that used in Ref. [23| but the expressions for the ma- 
trix elements have been improved. The basis functions 
arc localized piecewise polynomials defined on a prod- 
uct grid in the prolate spheroidal coordinates, with some 
matrix elements evaluated within the discrete variable 
representation (DVR) approximation [25rl28j. For odd 
A4 values the basis functions include (unitless) factors of 

p = ^/(£ 2 — — y l 2 ) to cn f° rcc square root boundary in which linear combination the second derivatives 
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FIG. 1: (Color online) Outgoing waves ip ac (u)) for hv=5.9, 14.6, 18.95, 23.3, 33.1, and 56.0eV, evaluated at r\ = 1, i.e., along the 
bond axis. The behavior of the wave function is trivial in the r\ coordinate. The modulus is plotted on a logarithmic scale, and 
the color denotes the phase, with the dark-medium-light (blue-red-yellow) colors having time derivatives light-dark-medium 
(yellow-blue-red). The coordinates are real valued within the plots. 



have been cancelled, and more trivially that 
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The raising and lowering operators are 
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As mentioned above, the primitive basis functions are 
defined with factors of p = y/ (£ 2 — 1)(1 — rf) for odd Ai . 
The matrix elements of the raising and lowering operators 
for a bracket with the ket having even quantum number 
M. therefore involve integrals of the following operator 
with respect to the polynomial basis functions 



{M ± i)£n 

(7) 



The matrix elements for odd-A'l ket are in turn expressed 
in terms of matrix elements of the operator 



(e - rp)l±p = ± 



(8) 



The matrix elements of the individual terms in Eq. @ , 
which are hermitian, of the operator in Eq. (|5j) and Eq.® 
that occurs both in rj and £, which is antihermitian, and 



of the antihermitian ( ^ — ) operator are integrated 
exactly by quadrature. As in Refs. [U, [23| only one di- 
mensional integrals are involved and therefore the matrix 
representations of these operators are quite sparse. 

The basis set employs exterior complex scaling j30l - [36j 
of the electronic and nuclear coordinates in order to en- 
force outgoing wave boundary conditions exactly. The 
coordinates of electrons and nuclei are rotated into the 
complex coordinate plane in the asymptotic region, which 
results in an antihermitian component of the Hamilto- 
nian that only absorbs outgoing flux. Bound states are 
not absorbed and despite the fact that Rydberg states 
penetrate into the complex scaled region, their analytic 
continuations are accurately represented, obeying an or- 
thogonality relationship, having unperturbed real eigen- 
values, etc. 



II. TOTAL CROSS SECTIONS 



The absorption cross section is calculated |36h38| a la 
Lewis-Dalgarno [39| by solving the linear equation 



(9) 



in which Eq is the energy of the initial ground rovibronic 
eigenstate oj is the photon energy, p is the dipole 
operator, and G + (Eq + to) is the outgoing wave Green's 
function as represented by exterior complex scaling. 

Examples of the calculated time independent half- 
scattering wave functions ip sc are shown in Fig. [TJ The 
wave functions are evaluated parallel or perpendicular to 
the bond axis for the corresponding laser polarizations, 



3 



10" 



10" 



10" 



CO 1 Q 3 



10" 




Optical total 
J?^ D.I. ♦ 




// D.E. ♦ 




^jf^ssssa^ggl 








A 











10 



20 



30 
hv (eV) 



40 



50 



60 



FIG. 2: (Color online) Cross sections calculated for total <j(oj), 
light colored lines, dissociative ionization a DI (ui), black small 
diamonds, and dissociative excitation a DE ((jj), large orange 
(pale) diamonds, as calculated via Eqs. (|10I16|I . and total via 
the optical theorem, Eq. (|17jl . thick black lines, for parallel 
(top) and perpendicular (bottom) polarization. 



i.e., at the point 77 = ±1 or 0, correspondingly. The 
behavior in the r\ is mostly uninteresting, being mostly 
p-wave outgoing flux for dissociative ionization, for in- 
stance. 

To extract the cross sections from the ip sc the method 
of Ref. [4(| as adapted to exterior complex scaling in 
Ref. [23[ is applied. As the outgoing wave at a given 
photon energy is directly calculated via Eq. |H1 no Fourier 
transform is needed. The total breakup cross section is 
obtained in the length gauge and for polarization parallel 
to the bond axis via 



<j{uj) = ^Traujh(4> sc (uj)\&(H)\i' sc {uj) 



(10) 



with a the fine structure constant. In this expression 
a(if) = ^(H—W) is the antihermitian part of the Hamil- 
tonian, the hermitian part being h(H) = ^(H + H>). 
This is an isotropic cross section so there is the factor 
of |. For perpendicular polarization there is a factor of 
I and so the corresponding coefficient in Eq. (fTU|) is -y. 
The result is shown in Fig. [2] 

To distinguish dissociative excitation from dissociative 
ionization the antihermitian part of the hamiltonian is 
divided into the part that absorbs flux for large bond 
lengths R and that which does so for large values of the 
prolate spheroidal coordinate £. With the identity 
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and the shorthand 
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the full Hamiltonian may be abbreviated 



V + DY 



(12) 



(13) 



and the antihermitian part of the Hamiltonian divided 
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such that 
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(15) 

The cross sections for dissociative ionization a DI (E) and 
dissociative excitation a DE (E) are calculated as 
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with a(u) = a DI (u) + a DE {w). 

The cross sections as thereby separated arc shown in 
Fig. [2j The gross features are apparent including peaks 
corresponding to the 1 and 2 H u or excitations. The 
total cross section is also calculated via the optical theo- 
rem, i.e. 



(j(uj) = ——iraujh im 



G+{E +lo) 



(17) 

The agreement between the two formally equivalent re- 
sults is exact although they are calculated quite differ- 
ently. 

The terms that occur in Eq. (|15[) are separately shown 
in Fig. [3] One can see that the terms are both posi- 
tive and negative. For both polarizations but especially 
for parallel, there is a negative contribution to the disso- 
ciative ionization cross section from the term h(D)a(Y) 
that cancels a part of that from the electronic hamilto- 
nian h(l/R 2 ).... The cross section at low energy that 
is entirely dissociative excitation is vastly dominated by 
the contribution from a(T), the antihermitian part of an 
operator that is very close to the hamiltonian for nuclei 
only. 

It is true that there may be outgoing flux that is ab- 
sorbed in the region in which both H^ nti and H^ nti are 
nonzero, calling into question the separation described 
above. However, the results for a DI and a DE presented 
here are converged with respect to the complex scaling 
radii in both the electronic and nuclear coordinate. 
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FIG. 3: (Color online) Contributions to the calculated cross 
sections from the antihermitian components of the complex 
scaled Hamiltonian as in Eqs. (|12ll3U15[) . The individual con- 
tributions are labeled in the figure and the total breakup cross 
section u(uj) is plotted with a bold black line. Top, parallel 
polarization; bottom, perpendicular. 
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FIG. 4: (Color online) Partial cross sections for dissociative 
excitation into final states of Rydberg H + H + . Cross sections 
a^ E (uj) are plotted with dotted lines and the sums of off- 
diagonal terms <x„ (ui),i 7^ j of Eq. (|18l) are plotted with dots 
only. Top, parallel polarization; bottom, perpendicular. 



III. DISSOCIATIVE EXCITATION 

The dissociative excitation cross section a DE is divided 
into contributions of the final electronic states of the hy- 
drogen atom. The division of the cross section proceeds 
via 

og E ( U ) = l™uh{r c {u)\uR)){UR)\* (18) 

such that = a DE (uj), where is the 

fixed nuclei electronic state. Both diagonal (i — j) and 
off-diagonal (i ^ j) contributions are calculated. The 
off-diagonal contributions are not expected to go to zero, 
because the prolate spheroidal coordinate R is not ex- 
actly the same as the dissociative coordinate, due to the 
mass of the electron, and this effect may not be reme- 
died by projecting at larger bond lengths. Presently the 
final states have been calculated in the Born Oppen- 
heimer approximation and also by using the improved 
adiabatic Hamiltonian of Ref. [|| without the coriolis cou- 
pling term, and the results are almost indistinguishable. 
In any case, the results should only be regarded as reliable 
when the diagonal cross sections, which should represent 
the actual cross sections, are much greater in magnitude 
than the off diagonal cross sections, which are presum- 
ably an artifact of the mixing of the two states involved 
by nonadiabatic coupling. 

Results are shown in Fig. 2J This figure shows the total 
cross sections for dissociative photoexcitation, summed 



over final states with a given principal quantum num- 
ber. In order to ensure an accurate assignment of the fi- 
nal electronic states, the projection upon the odd-parity 
states correlating to a given quantum number should be 
performed at a bond length large enough such that no 
crossings between these states and those correlating to 
other principal quantum numbers remain to be avoided 
as the bond length further increases toward infinity. In 
parallel polarization, the last avoided crossing involving 
an n = 3 electronic state occurs at approximately 16.5ao, 
and in perpendicular polarization the last such crossing 
for n = 4 occurs at approximately 15ao. The calcula- 
tions shown in Fig.|¥]thcrcfore have been performed using 
exterior complex coordinate scaling in which the scaled 
region begins at bond lengths greater than these. The 
results arc indeed not converged for smaller scaling bond 
lengths. These calculations are not large, but the com- 
puter code constructed to perform these calculations is 
deficient and therefore a bond radius sufficient to calcu- 
late, for instance, the n = 4 cross sections in parallel 
polarization has not been attempted. The last avoided 
crossing involving the n = 4 states occurs at approxi- 
mately 40ao. 

One may see that both the off-diagonal cross sections 
a ij E ' * 7^ i an d diagonal cross sections, i = j , both 
exhibit features at low energy that are congruent to the 
major peaks at low energy corresponding to the 1 E u 
and 1 Tl u states. Furthermore the magnitudes of the 
off-diagonal (unphysical) and diagonal (physical) cross 
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sections are comparable in this region, and the diago- 
nal cross sections are nonzero below the corresponding 
thresholds. These features therefore should be regarded 
as spurious, and probably relate to the fact that the 
asymptotic states with which the projection is performed 
are not precisely the true asymptotic states, as explained 
above. These results at low energy, those for n . = 2 and 
3 for parallel polarization and n = 3 and 4 for perpen- 
dicular polarization, are therefore spurious. 

At higher energy, there are significant off-diagonal (un- 
physical) cross sections for perpendicular polarization, 
but the magnitude of these cross sections are several or- 
ders of magnitude below the physical ones. The diagonal 
cross sections should be regarded as reliable in the high 
energy region. 



IV. DISSOCIATIVE IONIZATION 

The dissociative ionization flux at a given photon 
energy may be differentiated with respect to the en- 
ergy sharing between the electronic and nuclear degrees 
of freedom, and with respect to the angular behavior. 
Presently the angular behavior is not analyzed, and was 
found to consist of essentially p-wave correlation between 
the electronic and nuclear degrees of freedom for nearly 
all photon energies, as would be expected. The disso- 
ciative ionization cross section is obtained by calculating 
the amplitudes Ai(k, k) for breakup as a function of ki- 
netic energy sharing and electronic angular momentum 
quantum number, 



Mk,K) = (^-{k,K)\^\^°) 



(19) 



Presently amplitudes Ai(k,n) are calculated exactly, 
and also using two degrees of approximation, as described 
below. 



A. Exact and approximate amplitude expressions 

For the fixed-nuclei problem, exact final states were 
calculated in Ref. [2l[ by solving the equation 

*,-(£) = 4>o + G+ {E){H - E)^ (20) 

with the zeroth order wave function a coulomb wave </>o = 
fi(kr)Pi m (cos9), E = 4-, and the interaction term 



(H - E)(f,o 
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(21) 



In contrast, for three body scattering with pairwisc 
interactions Eq. (|20| is not valid. As an alternative to 
employing an explicit representation of ^~(E), station- 
ary phase expressions [U |42| that can be implemented 
in a numerically robust way [36l438| have been applied 
to three-body coulomb breakup problems with two elec- 
trons. 



However, the prolate spheroidal coordinate system, 
along with the unequal masses between the electronic and 
nuclear degrees of freedom, in general appears to allow 
Eq. (f2"0)) to be implemented such that it yields an accurate 
final state. At the end of the electronic grid in prolate 
spheroidal coordinates, the electron is always at a greater 
radius than the nuclei. Thus, we should not expect to be 
able to construct final states for which \x^k < fi e n, for 
which the protons recoil faster than the electrons and 
are thereby shielded from one another by the electron. 
Given a maximum of approximately 13.6eV nuclear ki- 
netic energy release, this would indicate that our results 
are certainly good above 7.5meV electron energy, a quan- 
tity that is not visible on any of the figures. 

I therefore solve 

¥,-(fc, k) = f t {kr)fi(KR) + G+(E)(H - E)f l {kr)f x {KR) 

(22) 

wherein fi(kr) and fi(nR) are attractive (Z = 2) and 
repulsive (Z — 1) Coulomb functions in the electronic 
and nuclear degrees of freedom, energy-normalized. This 
is performed in the same manner as in Ref. [21( . The final 
state wave functions so constructed are orthogonal to the 
bound rovibronic states for M. = and 1, to within no 
more than one part in 10 3 , in general. 

With the amplitudes defined as per Eq. ([T9| , cross sec- 
tions differential in the kinetic energy sharing between 



the electron 



e and that of the nuclei 



constant total energy E + e, are 



d_ 
dE 
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-it aurni 



E+e 



^2\Mk,K)\ 



E, at 



(23) 



j4,(M) = (tfj-(M)|0|¥°) 



wherein \& — {k, k) is energy normalized. 

Approximate final states are often employed in cal- 
culations in the literature, and in some contexts simple 
unperturbed product wave functions 0o are surprisingly 
accurate. For instance, in time-dependent calculations on 
small atoms and diatomics, cross sections may be calcu- 
lated [43Tj46j by projecting a propagated wave packet 
onto unperturbed coulomb wave functions, as long as 
enough time has elapsed such that the ionized electrons 
have escaped beyond the molecule. In systems contain- 
ing resonances, this method becomes less tractable the 
longer-lived the resonances are. A comparison of different 
amplitude expressions for single and double ionization of 
two electron systems, similar to that presented here for 
dissociative ionization of H^, can be found in Ref |47j . 

Approximate final states (k, k) are constructed as 
products 



*,-(*.«) ~ fi^R)^(k 2 /2-R) 



(24) 



of coulomb waves in the bond distance and the exact 
fixed- nuclei scattering state tp~ . This is therefore a Born- 
Oppcnhcimer representation of the final state. Finally, 
the initial state is replaced with its Born-Oppenheimer 
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FIG. 5: (Color online) Total cross sections for dissociative ion- 
ization calculated in different ways: Eqs. |16] [23] [25]) . PSIMI- 
NUS with that for dissociative excitation for reference. Top, 
parallel polarization; bottom, perpendicular. 

approximation as well, such that the amplitudes are the 
matrix elements of the Born-Oppcnhcimcr amplitudes 
with respect to the initial and final vibrational states, 
which expression is called the Chase approximation |48| : 

At(k,n) « ( dRh{KR) X o{R)Mk;R) , (25) 



in which \o is the ground vibrational Born Oppcnheimcr 
state. 



B. Results 

Cross sections for dissociative ionization arc shown in 
Fig. [5] For these total cross sections, there is very little 
difference among the various results calculated using the 
different amplitude expressions, exact and approximate. 
However, in perpendicular polarization, there is a large 
discrepancy between these results and <j di (lu) as defined 
by Eq ll6l The origin of this discrepancy is unclear, but it 
calls into question the division of the cross section as de- 
fined by that equation. Further study of this discrepancy 
is therefore indicated. I regard the results calculated via 
the amplitude expressions to be reliable, due to the fact 
that one of them has been calculated in a formally exact 
manner. 

The distributions of kinetic energy between the elec- 
trons and nuclei are shown in Fig. [6] In these figures the 
exact, approximate, and chase approximation results are 
compared. One can see that for the majority of the plot- 
ted area the results are in agreement, but at threshold 
there are significant disagreements. The chase approx- 
imation yields qualitatively incorrect behavior, yielding 
strong minima in the cross section as a function of energy 
sharing whereas in the exact result there are none. 
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